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Abstract 

Recent experiments showed that Co undergoes a phase transition from ferromagnetic hep phase 
to non-magnetic fee one around 100 GPa. Sinee the transition is of first order, a eertain region of 
co-existence of the two phases is present. By means of ab initio calculations, we found that the 
hep phase itself undergoes a series of electronic topological transitions (ETTs), which affects both 
elastic and magnetic properties of the material. Most importantly, we propose that the sequence 
of ETTs lead to the stabilisation of a non-collinear spin arrangement in highly compressed hep 
Co. Details of this non-collinear magnetic state and the interatomic exchange parameters that are 
connected to it, are presented here. 
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By the development of the diamond anvil cells, it became possible to carry out high pres¬ 
sure experiments upon to a few hundreds GPa.[l] Such compressions roughly correspond to 
the state of the matter close to the Earth’s core, so that the processes occurring under the 
deep mantle can be reconstructed in nowadays laboratories. [2, 3] Particularly, the investi¬ 
gation of magnetic systems is of interest, because it can give fundamental insights on the 
origins of the geomagnetism. [4] 

Cobalt is not a particularly abundant element in the Earth’s core, and is therefore not 
intensively studied under extreme conditions. However, it is an outstanding element from 
the technological and scientific viewpoints. At ambient conditions, Co crystallizes in an hep 
lattice showing a ferromagnetic (EM) order with a high Curie temperature (Tc) of 1388 K. 
Its Eermi surface (ES) is highly spin-polarised, which makes it a great spin-filter material. 
The hep phase can be stable up to a very high pressure, but at around 100 CPa an hep 
to fee transition takes place. [5, 6] There is an on-going debate about the magnetic state of 
cobalt above this critical pressure. Iota et al. proposed that Co remains in the magnetic 
hep phase at least up to 100 CPa. Eor higher pressures, hep Co is gradually transformed 
into the non-magnetic (NM) fee phase, which identifies a region of 50 CPa where the two 
phases coexist and Co gradually loses its magnetism. [7] On the other hand, Torchio et al. 
argued that the magnetism is already completely lost at 120 GPa in the mixed phase. [8] 
This is also in qualitative agreement with more recent results reported by Ishimatsu et al, 
but these authors also suggest the existence of a super paramagnetic (PM) fee phase in 
the high-pressure region above 135 GPa. [6] All aforementioned groups analysed the x-ray 
magnetic circular dichroism (XMCD) data, measured at the K edge. It is worth mentioning 
that this type of experiments does not probe directly the spin moment, but rather the weak 
moment of the 4p states. [9] Thus, some small values of the total magnetisation might appear 
as a negligible noise on the spectra. 

Theoretically, the problem of magnetism of Co under pressure has been extensively inves¬ 
tigated by means of the density functional theory (DPT).[10-15] Nevertheless, the definition 
of the transition pressure is still uncertain, which can be traced back to differences among 
the employed DPT implementations. [16] The most recent calculations indicate that the fee 
phase favours a NM state before the transition, while the hep phase remains magnetic up 
to more than 160 GPa. [17] The FM-NM transition in the hep phase was identified to be of 
second order, while in fee phase it is a first-order one. [12] 
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On the experimental side, the pressure-induced anomalies of the c/a ratio[18] and Raman 
frequencies[19] were attributed to a strong magneto-elastic coupling, which is argued to be an 
intrinsic feature of the material. However, any unequivocal evidence for the coupling between 
magnons and phonons have not been presented yet. In this Letter, we point out that the 
above-mentioned peculiarities result from an electronic topological transition (ETT), or the 
so-called Lifshitz 2| transition. [20] We also argue that the change of the FS gives rise to 
a non-collinear magnetic ground state in hep Co under pressure, and that this may be an 
excellent way to experimentally detect the influence of the ETT. 

In the present study we have used several DFT implementations. The structural opti¬ 
misation for both hep and fee phases was performed using the Vienna Ab initio Simulation 
Package (VASP).[21, 22] For this purpose, the stress tensor was analysed to reach a certain 
value of external pressure. The lattice parameters were converged by using a large plane- 
wave energy cut-off of 600 eV and very dense k-point grids, i.e. 45x45x27 and 45x45x45 
for respectively hep and fee structures. Then, the obtained crystal structures were used to 
perform DFT calculations by means of the full-potential linear muffin-tin orbital method 
(FP-LMTO) as implemented in the RSPt code[23]. These simulations were used to extract 
the effective exchange parameters (Jjj) and calculate the corresponding adiabatic magnon 
spectra, as well as for analysing the FS. The J^j’s were computed between the 3d states, 
which were projected onto the muffin-tin spheres (for details see Ref. 24). Additional sim¬ 
ulations of the spin spirals (SS) were performed via the PY-LMTO code. [25] All results were 
obtained with the GGA-PBE[26] exchange correlation functional. Within this framework, 
the experimental hep-fee transition pressure was shown to be well reproduced. [12, 13] We 
did not consider the effects of strong correlations, since it would require to know not only the 
value of the Hubbard U parameter, but also its pressure dependence. Moreover, in our study 
we have to deal with very high pressures, where the large band broadening should diminish 
the correlation effects. Spin-orbit coupling is neglected throughout the paper, however we 
discuss how it possibly affects the results. 

The calculated FS of the ferromagnetic hep Co as a function of pressure is shown in Fig. 1. 
At ambient conditions, four minority bands and two majority-spin ones cross the Fermi level 
{Ep). The shape of the FS is barely changed up until roughly 80 GPa (see Supplemental 
Material (SM)[27] for more details). Above 80 GPa, a sequence of ETTs happen at different 
pressures in the spin-up channel: at around 80, 100 and 140 GPa, one more band is counted 
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to contribute to the Ep. At 180 GPa the NM state is reached, which is characterised by all 
five bands in each spin channel crossing Ep. 



FIG. 1. (Color on-line) Fermi surface of the majority-spin (f) and minority-spin (j,) electrons in 
hep Co for different values of external pressure (in GPa). The results were obtained with the 
GGA-PBE functional. The plot was produced with the XCRYSDEN[28] software. 

Close to 80 GPa a particularly important ETT takes place. Many small sheets emerge in 
the spin-up channel, and the part of the FS close to the P point of the spin-down channel is 
also significantly modified. This causes a transition from a strong ferromagnet to a weak one. 
Although there is no direct experimental evidence of this ETT, the changes in the topology of 
the FS can explain the changes that have been observed in several other quantities at around 
80 GPa. First of all, at around this pressure, the cja ratio alters its pressure derivative, as was 
unveiled by Antonangeli et a/, [18] who confirmed it both experimentally and theoretically. 
Second, vibrational properties, such as e.g. sound velocities, show anomalous behaviour at 
around 75 GPa. [29] This anomaly is accompanied by a change of the slope of the E 2 g phonon 
modes, which was first measured by Goncharov et al. (Ref. 19) and later on reproduced in 
the DFT calculations.[13] Moreover, previous DFT studies[14, 17, 18] suggested that at the 
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same compression ratio, i.e. VjV^ ~0.8, the pressure-driven reduction of the spin moment 
(Ms) becomes faster. These anomalies in the c/a ratio and Ms are very well reproduced 
in our calculations, as reported in the SM.[27] Most importantly, in combination with our 
results on the FS, they can offer a strong evidence for the existence of the ETT, and can be 
seen as signatures of this order transition. At low pressure, hep Co has a majority-spin 
band completely occupied, which results in a weak pressure dependence of Ms- However, the 
ETT pushes up the majority-spin band towards Ep, so that Ms becomes more vulnerable to 
external stimuli, e.g. to pressure (for a detailed explanation of this mechanism, see Ref. 30). 
The resulting decrease of the Ms is the primary reason for the destabilisation of the EM 
hep phase at high pressures. This is similar to the bcc-hcp transition in Fe, which can be 
explained by mapping of the DFT solution on a simple Stoner model. [31] Once the energy 
gain due to magnetic moment formation, which is proportional to (Ms)^, gets too small, the 
system transforms into more close-packed structure. 



FIG. 2. (Color on-line) Calculated elastic constants (upper panel) along with the DOS at the Ep 
(lower panel) as a function of pressure. DOS is summed up for both spin channels. 

Signatures of the ETT can also be found in the elastic constants, which are reported in 
Eig. 2 together with the density of states (DOS) at the Ep. First of all, one can see that 
below 70 GPa DOS(£'^) is weakly dependent on the applied pressure. At 80 GPa, a small 
hump is visible, which is a manifestation of the ETT, as discussed previously. After the ETT, 
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the DOS{Ef) becomes strongly pressure-dependent, and quickly surges upon compression. 
As illustrated in Ref. 32, there is clear link between ETT and lattice properties, and the 
elastic constants result to be directly related to the derivative of the DOS at the Ep. In the 
case of hep Co we observe that all elastic constants except 044 change their trends around 
80 GPa. Similar findings have already been reported in a previous study by Kuang et a/. [15] 
But from our results it becomes apparent that this happens as a result of a remarkable 
change of the DOS(£'i 7 '), associated with the ETT, which appears at the same pressure. 

In the literature, these anomalies in the magnetic, elastic and vibrational properties of 
Co were associated with a magneto-elastic coupling, without any clear explanation. [13, 14, 
18, 19, 29] In our work we provide a strong evidence that all these changes are consequences 
of the ETT. Present scenario is analogous to the one proposed for PM hep Fe.[33] In Ref. 33 
it was reported that only way to observe the ETT was to augment the DFT solution with 
the dynamical mean field theory (DMFT)[34]. DFT-|-DMFT[35] significantly improves the 
description of the PM state as compared to the bare DFT, due to its proper account of the 
local moment fluctuations. The latter ones give rise to a strong damping of the quasiparticles 
and the pronounced correlation effects are expected. Our study, however, concerns the 
ordered phase where the quasiparticles close to the Ep have much longer lifetimes and the 
many-body effects are moderate. [36] 

Note that ETTs do not necessary lead to the anomalies, for instance, in the c/a ratio. 
An example of such situation is hep Os. One of the reasons for this is that the value of 
the DOS{Ep) is smoothly dependent on the pressure even when the ETT occurs.[37] In 
our case these changes are more conspicuous (see Fig. 2 ). Another reason is probably a 
weak coupling between the electrons and the lattice. On the contrary, it was recently shown 
that the electron-phonon coupling in hep Co is five times larger than that in bcc Fe.[38] In 
our opinion, a combination of these two factors in hep Co allows the ETT to influence the 
lattice-related properties of the material. 

In the range of 80-180 GPa the hep phase of Co remains magnetic and its spin moment 
strongly depends on the volume (see SM[27]). However, the stability of the EM state at such 
conditions has never been examined. To shed light on this, we have calculated the effective 
Jjj-parameters by means of the magnetic force theorem. [39] The exchange parameters within 
the first few neighbouring shells, calculated at various pressures, are shown in Fig. 3 (upper 
panel). 
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FIG. 3. (Color on-line) Upper panel: Inter-site exchange parameters in hep Co as a function of 
distance for several values of pressure. Lower panel: Calculated adiabatic magnon spectra, obtained 
using the Jjj’s shown in the upper panel. The Jjj’s are extracted for the Heisenberg model with 
unitary spins (see e.g. Ref. 24.) 

The zero-pressure values are in good agreement with prior DFT studies. [40, 41] When 
going from 0 to 80 GPa, one can see that the nearest-neighbour couplings (Ji) remain rather 
close to each other. This difference is instead most pronounced for the Jq interaction: it 
is antiferromagnetic at ambient conditions and gets strongly enhanced at 80 GPa. This is 
the primary reason for the emergence of magnon softening along the P — A direction, as 
illustrated in the adiabatic magnon spectrum reported in the lower panel of Fig. 3. Further 
compression leads to vast modifications of the exchange parameters. Ferromagnetic Ji and 
J 3 couplings become suppressed and the latter one eventually changes sign. At 140 GPa the 
Mg is about 0.7 hb per atom, i.e. the system is still magnetic, but its FM state is clearly 
unstable. The global minimum is formed at around q = (0,0,0.237r/c), which is along F — A 
(or z) direction. 
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The instabilities in the magnon spectra, obtained under a certain pressure, indicate that 
the long-wavelength excitations destroy the FM ground state. What is peculiar is that they 
have only been found along the F — A path in the Brillouin zone (BZ), which motivated us 
to simulate the single-g family of the magnetic states. Hence, we have calculated the total 
energies of a set of the SS conhgurations by means of the generalised Bloch theorem. [42] 
Any SS state is characterised by two parameters: propagation vector and the cone angle 
between the direction of the propagation and the magnetisation (0). In Fig. 4 we show the 
relative energies of the SS states with respect to the FM one. 



FIG. 4. (Color on-line) Calculated total energies (relative to that of the FM state) of the spin 
spiral states propagating along F-A direction. At each pressure the calculations for two different 
cone angles (©) were performed. Dashed lines correspond to ©=30 degrees, solid ones to ©=90 
degrees. No peculiarities were found in other directions and thus they are not shown. 

An inspection of Fig. 4 suggests that at low compressions (e.g. 30 GPa) any SS state 
has a higher energy than that of the FM solution. This is consistent with the stable spin 
wave excitation spectra in this regime. However, this behaviour is changed as the pressure is 
raised up to 90 GPa. In a certain range of q vectors the SS states become energetically more 
favourable. The family of states with 0=90 degrees seem to be more stable than all other 
considered states. Further increase of the pressure shifts the total energy minimum towards 
higher g-values. For instance, at 130 GPa the minimum becomes even more pronounced 






and the total energy difference with respect to the FM state reaches 0.05 mRy. We have 
also performed the calculations for the spirals propagating along other directions in the BZ, 
but no solutions having lower total energy were found. Hence the 2 ; direction is particularly 
favourable for the propagation of the SS’s. The reasons for it can be seen from the analysis 
of the band structure, presented in SM.[27] 

Thus we suggest that the non-collinear states should be present in hep Co upon com¬ 
pression. In fact, it was previously proposed that it should be possible to stabilise the SS 
configurations in all ferromagnetic transition metals upon compression. [43] However, this 
situation is very seldom observed experimentally, due to the presence of the structural tran¬ 
sitions and/or the collapse of the magnetic moment at lower pressures. According to Ref. 43, 
the occurrence of the non-collinear spin states is governed solely by the FS topology. In par¬ 
ticular, it is absolutely necessary to have both spin-up and spin-down bands crossing Ep. 
We believe that hep Co under pressure is a prototype system, where the ideas, proposed in 
Ref. 43, are realised. As we have shown above, at around 80 GPa one of the spin-up bands 
start to cross the Fermi level, which leads to the stabilisation of the SS state. It was recently 
shown that under certain pressurizing conditions the pure hep phase of Co can be preserved 
up to at least 160 GPa. [44] Thus the proposed non-collinear states should exist in a wide 
pressure range, which in principle leaves enough room for their experimental detection. 

Our results indicate that the behaviour of cobalt under pressure is in many ways remi¬ 
niscent to that of hep Fe. In the PM state of Fe the c/a anomaly was shown to originate 
from the pressure-driven ETT.[33] Moreover, the emergence of non-collinear SS states was 
also suggested from the analysis of the band structure. [45, 46] Extensive investigation of 
hep Fe was partially motivated by the discovery of superconductivity in this system below 
1.5 K.[47] A close similarity of the two systems implies that an experimental investigation 
of the low-temperature properties of hep Go would be of high interest. 

In relation to the K-edge XMCD experiments[6, 8] it is worth mentioning that a 90-degree 
SS state does not have any magnetic dichroism, because all components of the magnetisation 
are compensated on a macroscopic scale. On the other hand, an application of the external 
magnetic field is necessary in these experiments. Large fields, in principle, could lead to the 
bending of the SS, thus resulting in a different cone angle of the spiral. Our results suggest 
that, even if 0 reaches 30 degrees, these SS states will be more favoured than the FM one 
(see Fig. 4). However, such a state would already be characterised by a non-zero XMGD 
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signal, since the total Mg will be finite. 


At this point it is important to discuss the relativistic effects, so far neglected in our study. 
Most importantly, spin orbit coupling gives rise to the magneto-crystalline anisotropy energy 
(MAE), which defines a preferable direction (or plane) for spins to point to. At ambient 
conditions the MAE of hep Co is 4.77 /rRy per atom. [48] This value is one order of magnitude 
smaller than the maximal difference between the total energy of the SS and the EM state. 
Moreover, the most favourable SS configuration was found above 100 GPa. At such large 
compressions the orbital moment is strongly suppressed, since it attenuates faster than the 
spin moment. [9, 44] The MAE is supposed to follow the orbital moment, since their relation 
is known. [49] Thus, we are confident that our conclusions remain valid in the case of finite 
spin-orbit coupling. 

In summary, we have shown that a sequence of ETTs happen in hep Co upon compression. 
The most important transition, occurring at 80 GPa, is suggested to be the reason for the 
number of anomalies in the magnetic, elastic and vibrational properties, confirmed both 
experimentally and theoretically. In the high-pressure regime, where the magnetic moment 
of Co is reduced to ~ l.S/i^, the ETT is demonstrated to lead to the stabilisation of non- 
collinear magnetic states. 
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SUPPLEMENTAL MATERIAL 


Band structure under pressure 

The Fermi surface (FS) changes, presented in the main text of the article, clearly show 
the presence of the series of electronic topological transitions (ETTs) in hep Co upon com¬ 
pression. Another way of observation of the ETT is an analysis of the band structure. In 
Eig. 1 we show the fc-resolved spectral function for different pressures. 

An inspection of Eig. 1 reveals that the most of the ETT’s occur for spin-up bands. Above 
80 GPa there are several transitions particularly along M-L direction in the Brillouin Zone, 
which is parallel to the .2 axis. These transitions create a larger overlap between majority 
and minority states. As a result, a particular family of spin spirals, which is characterised 
by q parallel to z, becomes preferable at high pressure. 

ETT and other properties 

In this section we provide a direct evidence that the ETT’s in compressed hep Co give 
rise to the anomalies in c/a ratio and magnetic moment (Ms). In Eig. 2 the pressure depen¬ 
dence of both quantities are presented along with the corresponding values of the density of 
states (DOS) at the Eermi level. The latter quantity varies highly non-monotonously under 
pressure, because of the sequence of the ETT’s happening upon compression. However, the 
most remarkable changes in the DOS(£'^) and its pressure derivative occur around 80 GPa. 
This is the main reason why this particular ETT affects so drastically the properties of the 
material. 

The trends of the magnetic moment and c/a ratio are very similar to the ones obtained 
by Antonangeli et al. (Ref. 18). The differences in the reported results are more pronounced 
for the c/a ratio and clearly originate from the different basis sets used in the calculations. 
The authors of Ref. 18 employed full-potential linearized augmented plane wave method 
(FP-LAPW), while we used projector augmented wave method (PAW). Moreover, c/a ratio 
is a parameter, which is hard to converge. Its final values might depend on the method for 
its extraction. We used the stress tensor calculations instead of the parabolic fit of total 
energy as was done in Ref. 18. It is also worth emphasising that we have used a much denser 
fc-point grid containing 45x45x27 points, as compared with Antonangeli et al. who used 
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FIG. 5. (Color on-line) Computed fc-resolved spectral functions of the ferromagnetic hep Co at 
different pressures. Spin-up bands are shown on the left side of the plot, while spin-down states are 
on the right one. The results were obtained with the RSPt code. Red circles underline the bands, 
which cross the Fermi surface at different compressions. The Fermi level is set to zero. 

16 X 16x8 points. 

In spite of all these details, both reported trends are qualitatively the same and show 
anomalies around 75-80 GPa. Fig. 2 re-considers these results in the view of an ETT, which 
we observed. 
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FIG. 6. (Color on-line) Pressure dependence of several ground state properties. Upper panel: DOS 
at the Fermi level. Middle panel: optimised cja ratio. Lower panel: Magnetic moment {Ms) per 
Co atom. 


16 






